! ---- Choice of the time step  
 
subroutine dt_adjust
include 'precision.inc'
include 'params.inc'
include 'arrays.inc'


!--------------------------------------------------------------------------
!  G. Ito commented out this part of code because
!  we have set amul=1.  Also remesh uses min(dt_elastic,dt_maxwell)
!--------------------------------------------------------------------------
!dt_min = min(dt_elastic, dt_maxwell)	!G.Ito

! Adaptive scaling (Cundall, 1982)
!if ( boff .lt. ratl ) then		!G.Ito
    ! Increase the time step 
!    dt = dt * amul			!G.Ito
!elseif ( boff .gt. ratu ) then
    ! decrease time step
!    dt = dt / amul			!G.Ito
    !dt = max (dt, dt_min )
!endif      				!G.Ito

!write(*,'(a24,4E12.4)') 'elast, maxwell, min, dt=', dt_elastic, dt_maxwell, dt_min,dt  !debuggin

!dt = max (dt, dt_min )		!G.Ito

dt = min(dt_elastic, dt_maxwell)	!G.Ito

return
end
